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Abstract 

In treating the relativistic 3-quark problem, a dressed-quark propagator parame- 
terization is used which is compatible with recent lattice data and pion observables. 
Furthermore 2-quark correlations are modeled as a series of quark loops in the scalar 
and axialvector channel. The resulting reduced Faddeev equations are solved for 
nucleon and delta. Nucleon electromagnetic form factors are calculated in a fully 
covariant and gauge-invariant scheme. Whereas the proton electric form factor Ge 
and the nucleon magnetic moments are described correctly, the neutron electric 
form factor and the ratio Ge/Gm for the proton appear to be quenched. The influ- 
ence of vector mesons on the form factors is investigated which amounts to a 25 % 
modification of the electromagnetic proton radii within this framework. 
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1 Introduction 



In tackling the covariant bound state problem in QCD, models based on a combined Dyson- 
Schwinger (DS) and Bethe-Salpeter (BS) approach have found widespread application, for 
recent reviews see refs. [I], |2j. This approach has been most successful in describing light 
pseudoscalar mesons and their electromagnetic properties. Starting from a suitable model for 
the gluon and the gluon-quark vertex in the infrared, and using this model consistently for the 
q — q scattering kernel in the meson BS equation, these mesons retain their character as both 
q — q bound states and Goldstone bosons. Masses |§, decay constants [|J and form factors || 
are found to be in excellent agreement with experimental data. In these studies the so-called 
rainbow-ladder approximation is used which consists in retaining the bare quark-gluon vertex 
and a gluon propagator. The latter is modeled with an enhancement at intermediate momenta 
which provides enough strength to generate a dynamical quark mass. 



Along these lines the nucleon's bound state amplitude can be obtained by solving a relativistic 
Faddeev equation which needs as input the full solution for the q — q scattering kernel. It is 
known that in the rainbow — ladder approximation this kernel exhibits diquark poles |6], [/J, with 
scalar (0 + ) diquarks (~ 0.7 — 0.8 GeV) and axialvector (1 + ) diquarks (~ 0.9 GeV) having the 
lowest masses. Other diquark correlations have much larger masses. Although these poles might 
correspond to unphysical asymptotic states (and indeed disappear when going beyond rainbow- 
ladder 1, i HI) they give us a hint that + -l + quark-quark correlations are expected to be 



dominant in the nucleon. This argument receives support from recent lattice calculations [1 I 



and also explains the u—d valence quark asymmetry observed in deep inelastic scattering 12, 131. 



The full relativistic Faddeev problem is highly involved and has been solved so far only for 
a NJL model in lowest order where the q — q interaction is pointlike and therefore separable 
|T^| . If the q — q scattering (or t) matrix is separable, the Faddeev equations reduce to a 



quark-diquark BS equation which can be solved exactly. Inspired by this idea, the t matrix has 
been modeled in such a fashion in ref. [nj. Retaining free massive quarks and + -l + diquarks, 



electromagnetic, strong and weak nucleon form factors have been calculated, in good agreement 
with experiment (except for the magnetic form factors). It is noteworthy that especially the 
neutron electric form factor is positive and different from zero which is in contrast to the valence 
quark contributions in many non- or semi-relativistic quark models. (Due to the approximate 
SU(6) symmetry of these models, it is consistent with zero almost by construction.) Despite 
the positive results of the above mentioned models, the assumptions of free massive constituent 
quarks and diquarks is certainly too simplistic from a QCD point of view. Another line of 
approach has been taken by ref. ||16|| . In this study, a well-constrained parametrization of the 
quark propagator is used which was obtained by fitting it to a number of soft and spacelike 
meson observables [17|. It exhibits the basic feature of DS solutions: a mass function M(p 2 ) 
which is of the order of 400 MeV in the infrared and which evolves into the perturbative limit 
for p 2 — > oo. Furthermore this parametrization has no poles thereby mimicking confinement via 
the absence of a Lehmann representation. Scalar diquarks and the dominant nucleon Faddeev 
amplitude have also been modeled with entire functions (i.e. pole-free) and the electromagnetic 
form factors have been calculated. The results (fitted to Ge of the proton) show also a positive 
neutron Ge and enhanced magnetic moments due to the dressed quark propagator which also 
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leads, by use of the Ward-Takahashi (WT) identity, to a dressed quark-photon vertex. A 
drawback of this study is the lack of manifest electromagnetic gauge invariance. In order to 
maintain it the calculation of currents between bound states has to proceed by gauging (i.e. 
minimal coupling of the photon) to an interaction kernel and sandwiching the result between 
bound states which are solutions to bound state integral equations with exactly the same kernel 



18L 19 



In both studies [T5|, [L6| the parametrization of the q — qt matrix bears no relation to the quark 
propagator and thus to the dynamics which causes chiral symmetry breaking. There remains 
the possibility that the good results especially for the neutron Ge are rather a result of a clever 
parametrization of the q — q correlations than they reflect the underlying physics. Besides the 
bulk of contributions to observables coming from a quark core, one would also expect corrections 
to these mainly coming from the pion cloud. Their non-negligeable magnitude is apparent in 
recent lattice extrapolations to small quark (or pion) masses [f20| , |21|j , also recent covariant 



studies [|22| , [23|1 confirm that the nucleon mass shift due to pions is at least —200 MeV, thereby 



indicating the percentage level of pionic corrections to nucleon observables. 

Therefore, we will present an extension of the quark-diquark picture which, besides covariance 
and gauge-invariance, aims to include several constraints which are available through lattice and 
other QCD-phenomeno logical studies. Thus the number of free parameters will be confined, in 
fact to one, and we are in the position to explore the limits of a covariant nucleon quark core 
picture. We start from the main assumption to neglect three-quark irreducible interactions 
to arrive at solvable Faddeev equations. Evidence for this assumption is admittedly scarce, 



only in the limit of static quark sources lattice data [24]] seem to confirm a picture where 
flux-tubes between each pair of the three quarks minimize the free energy of the three-quark 
system. Proceeding from this assumption, we employ the above mentioned efficacious quark 
propagator parametrization which captures the essentials of the infrared behaviour of quarks 
within QCD to calculate separable + and 1 + diquark correlations by summing quark loop 
polarization diagrams (Sect. |2|). These correlations are employed to solve the nucleon and delta 
Faddeev equations (Sect. |]). Parameters are fixed by the masses of nucleon and delta, leaving 
only one free parameter which is essentially the extension of the diquarks. Form factors are 
calculated in a manifestly gauge invariance preserving scheme. Here, the construction of the 
diquark correlations ensures that the photon properly resolves the diquark. The dependence of 
the form factors on the diquark width is investigated (Sect. |J). Finally, in Sect. [S], we draw our 
conclusions. 

Throughout this paper we work in Euclidean metric (g^ u = 5^, {7^,7"} = 25^, 7 M ^ = 7 M ). 



2 Diquark correlations 
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Quark Mass Function 



Quark Renormalization Function 




Figure 1: The quark mass function (left panel) and the renormalization function compared to 
lattice data. 



2.1 The quark propagator 



Central ingredient to all calculations is the form of the quark propagator, 

S(p) = «7 -pa v (p 2 ) -a s (p 2 ) , (1) 



2\ 



ij>A{p 2 ) + Sip 2 ))' 1 . (2) 



if 1 + M(p 2 ) 

For the scalar and vector part we use the algebraic parametrizations: 

a s (x) = 2mJ(2(x + m 2 )) (3) 
+F{b x x) Hhx) [b + b 2 F{ex)] , 

a v (x) = [l-^(x + m 2 ))] , (4) 

with T{y) = (1 — e~ y )/y, x = p 2 /X 2 , rh = m/X, &s(x) = X&sip 2 ) and &v{x) = X 2 ay{p 2 )- The 
mass-scale is A = 0.566 GeV, and the parameter values are given by 



m 


bo 


h 


b 2 


h 


0.00897 


0.131 


2.90 


0.603 


0.185 



(5) 



In Fig. |T| we show the quark mass function M(p 2 ) and the renormalization function Z(p 2 ) for 
spacelike p in comparison with recent lattice data that have been obtained in Landau gauge 
|25| . Although the quark propagator fit has been performed to a number of meson observables 



within the DS framework [|17| and not to lattice data, the chosen parametrization represents the 
qualitative behavior of both functions very well. It is still too premature to ask for quantitative 
agreement since lattice calculations are not feasible for current quark masses around 10 MeV 
yet. We note that the lattice data indicate that the slope of the decreasing mass function is 
somewhat less steep than in the parametrization. We will find that this slope has influence on 
the ratio Ge/Gm of the proton. 
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Both functions a$ an d o"y are parametrized with entire functions. Thus they have no poles 
and reflect confinement. A major drawback, though, are the essential singularities at timelike 
infinity (p 2 = — oo). Consequently the quark renormalization function blows up for timelike 
momenta, and it has been shown in ref. [^| that this has disastrous consequences if one attempts 
to describe production processes where timelike momenta 0(1 GeV) are deposited onto the 
nucleon (Z(— lGeV 2 ) > 10 5 ). Nevertheless, for the bound state calculations described here, 
the quark propagator is needed at complex momenta where always \Z\ < 1. Due to technical 
obstacles we will calculate form factors only for Q 2 < 2 GeV 2 , and for these calculations the 
quark propagator is sampled at momentum points where \Z\ < 1.2. Thus for our calculations 
we do not expect artefacts of the parametrization to show up in the numerical results. 



2.2 The q — q t matrix 



According to the arguments given in the Introduction, we expect scalar and axialvector q — q 
correlations to be the most important ones within the nucleon. Thus we model a separable t 
matrix by 

t{k a ,kp] Pa ,pp)=t{k,p,P) = xld(k,P)D(P)x 5 jS (p,P) + (6) 

X ^(k,P)D^(P)r 7 s(p,P). (7) 

The relative momenta are defined as 

k M = ^(MpJ - MM) , (8) 

and the total diquark momentum is 

P = Pa + Pp = k a + kp . (9) 

We assume that the Dirac structure of the vertices x 5 (scalar diquark) and x M (axialvector 
diquark) is described by their leading components which non-relativistically correspond to 
quarks being in a relative s state: 

x%(p) = g 0+ (i 5 cUF(p 2 /w 0+ ), (10) 
x5s(p) = gMrcUF( P 2 /w 1+ ) . (ii) 

The scalar function T (defined below eq. (^)) describes the extension of the diquarks in their 
relative momentum variable, regulated by the widths Wq+ and W\+. Although the choice of this 
function is somewhat arbitrary, numerical results depend only on the diquark widths and not 
on the specific form chosen as we have checked by employing both monopole and dipole forms. 
The constants g$+ and g\+ are normalization constants yet to be determined. 

Antisymmetry between the quarks dictates the color and flavor quantum numbers of the vertices 
X- They are both in a color antitriplet representation. As (75C) is antisymmetric, the scalar 
vertex must also be antisymmetric in flavor space, in contrast to the axialvector vertex which 
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Figure 2: The separable t matrix and the definition of the diquark propagator 

is symmetric in flavor space due to the symmetric matrices (j^C). We adopt the following 
normalizations, restricting ourselves to the isospin subgroup: 

/_ \ P ABC 

{X» C )T = (12) 

„ \ C ABC 

f^H,C\AB _ { T 2n)ab e , . 

Capital letters denote color indices and small letters isopin indices. The r k represent the usual 
Pauli matrices. We will suppress color and flavor indices in the following, corresponding traces 
will have been worked out. 

We model the inverse diquark propagators D~ l and (D 11U )~ 1 by quark polarization diagrams 
as shown in Fig. |2| with an additional constant offset A + for the scalar diquark and A]+ for 
the axialvector diquark. This assumption describes the propagation of the quark pair being 
determined by an infinite series of loops as in Fig. 0. Similar expressions are obtained in the 



bosonized forms of the Global Color Model 27 or in the NJL model 14 . In the latter the 



constant A corresponds to the inverse strength of the four-quark interaction. 

For the scalar channel the inverse propagator reads 

D-\P 2 ) = -A^ + -n(P 2 ), (14) 
n(P 2 ) = -J ^Tr xV) S(P/2 + q) X 5 (q 2 ) S T (P/2 - q) . (15) 

The polarization function U(P 2 ) has also an essential singularity at timelike infinity (as the 
quark propagator). Typically it evolves from — oo at P 2 = —oo monotonically to zero at 
P 2 = oo, save for a tiny bump at around 4 GeV 2 . The effect of the constant A^ + is to shift D^ 1 
downwards that it acquires a zero at some — P 2 = m^ + . Thus, the propagator has a pole. This 
is very similar to the rainbow-ladder truncation of the quark DS eq./diquark BS eq. where 
poles do appear in the t matrix. At the pole we demand unit residue, 

A n(p 2 )| p2= _ m 2 =i, (16) 



cLP 2 v /l ^=- m o+ 

thus relating the hitherto unknown constant g + to m + or, equivalently, A +. 
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For the inverse propagator in the axialvector channel, we employ an ansatz similar to the scalar 
channel: 

(zr 1 r = -Al - IC ■ (17) 

The polarization loop 

n^(P 2 ) = - J ^Tr r (q 2 ) S(P/2 + q) X »(q 2 ) S T (P/2 - q) (18) 

can be split into longitudinal (pi*P u /P 2 n ax , L (P 2 )) and transverse ((5» u -P^P U /P 2 ) U^P 2 )) 
components. We note that IT aXj L(0) = IT aXj T(0) as it should be (to have no pole in the propagator 
at P 2 = 0). Furthermore we see from the numerical results that in the region P 2 e [—0.7, 1.5] 
GeV Il aXj L(-P) is approximately constant (deviations are ~ 1 %). Therefore we take the inverse 
propagator as 

(pnpv\ pvpv 
^-^r) -n ax , T (o)-^. (19) 

This is in accordance with the requirement that the longitudinal part of a spin-1 propagator 
not be dressed. 



The behavior of the transverse polarization, UaxT, is very similar to the scalar polarization, II, 
thus upon shifting by A 2 ^ the propagator acquires a pole at a mass which is larger than the 
scalar diquark mass if A + ~ A x +. At the pole a similar condition to eq. (|16D holds: 

1 n ax , T (p 2 )| p2= _ m 2 =i, (20) 



^p2 "ax.J-v /i±"=-m- + 

We wish to relate the constants A + and A x + which represent inverse coupling strengths in 
the scalar and axialvector channels, respectively. Consider a quark vector current-current 
interaction where the currents are color octet and Lorentz diagonal which arises e.g. from 
bosonizing a Global Color Model with the gluon propagator in Feynman gauge, ~ 5^ u . Upon a 
Fierz transformation into the scalar and axialvector diquark channels we find the relation |f2H| 



^ = 2 . (21) 
9i+ 9o+ 



In summary, we have parametrized the q — q correlations close to the rainbow-ladder truncation 
scheme which proved to be successful in the meson channels. Yet, scalar and axialvector 
diquarks are mainly characterized by an (unphysical) mass which should only be interpreted as 
an inverse effective correlation length than as a physical particle's mass. The composite nature 
of the diquarks is reflected by the propagators which are a series of quark loops. We used six 
parameters (diquark widths Wi, diquark normalization constants gi and quark-quark inverse 
coupling strengths Aj, i = {0 + , 1 + }) which are reduced to three free parameters by using the 



relations (p!6|,|20|,pTD . In the following, we will impose two more constraints using the masses of 



nucleon and delta, fixing essentially the inverse coupling strengths. 
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3 Faddeev equations for nucleon and delta 



A full derivation of the Faddeev equations for iV and A using separable q — q t matrices can 
be found in ref. [p9fl . In the following, we will only introduce the necessary elements. For 
the case of separable t matrices, it is convenient to introduce baryon-quark-diquark Faddeev 
amplitudes. These Faddeev amplitudes have to be decomposed in Dirac and Lorentz space 
after their projection onto positive energy states with spin 1/2 (JV) or spin 3/2 (A). 



3.1 Nucleon 



The nucleon Faddeev amplitude (or wave functions) can be described by an effective multi- 
spinor characterizing the scalar and axialvector correlations, 

*(p, PHP, s) = (fj*> P J] u(P, s). (22) 



V(p,P) 

u(P, s) is a positive-energy Dirac spinor (of spin s), p and P are the relative and total momenta 
of the quark-diquark pair, respectively. The vertex functions are defined by truncation of the 
legs, 

$m) = S ^ ( o (D^)- 1 ) (*") • (23) 

The coupled system of Faddeev equations for the nucleon wave and vertex functions can be 
written in the following compact form, 

/(2^ G " (P ' fc ' P) (^') (fc ' P)=0 ' ^ 

in which G~ 1 {p, k, P) is the inverse of the full quark-diquark 4-point function. It is the sum of 
the disconnected part and the interaction kernel. 

Here, the interaction kernel results from the reduction of the Faddeev equation for separable 
2-quark correlations. It describes the exchange of the quark with one of those in the diquark 
and thus the Faddeev equation reduces to an effective quark-diquark BS equation. Thus, 

G-\p,k,P) = (2vr) 4 5\p - k) S'\p q ) ^] Pd) ( ^)°-i (pd) )- 



1 / -xhvl) s T ( q ) xhvl) x"'(pi) s T ( q ) x 5 (pi) 

2 V x 5 (p 2 2 ) S T ( q ) S T (q) 



(25) 



Herein, the flavor and color factors have been taken into account explicitly, and x 5 ; X 11 stand for 
the Dirac structures of the diquark-quark vertices, see eqs. ( |IT||IT|) . The freedom to partition 
the total momentum between quark and diquark introduces the parameter rj G [0, 1] with 
p q = rjP + p and pd = (1 — rf)P — p. The momentum of the exchanged quark is then given by 
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q = —p — k + (1 — 2t])P. The relative momenta of the quarks in the diquark vertices x an d X 
are p 2 = p + k/2 — (1 — 3rj)P/2 and pi = p/2 + A; — (1 — 3i])P/2, respectively. Invariance under 
(4-dimensional) translations implies that for every solution \I/(p, P; 771) of the BS equation there 
exists a family of solutions of the form ^(p + (772 — f]i)P, P] 772)- 



Using the positive energy projector with nucleon bound-state mass M n , 

If 



A+ 



1 



the wave function can be decomposed into their most general Dirac structures, 



^(p,P) 



(S 1 + Jf JS 2 )A\ 



l"(A 3 + — i>A A ) lb k- 



(26) 

(27) 
(28) 



+ 



J?—(A 5 + — MebsA" 



In the rest frame of the nucleon, P = (0,iM n ), the unknown scalar functions Si and A± are 
functions of p 2 = p^p^ and of the angle variable z = P ■ p, the cosine of the (4-dimensional) 
azimuthal angle of p^. Certain linear combinations of these eight covariant components then 
lead to a full partial wave decomposition, see ref. |30| for more details and for examples of 
decomposed amplitudes assuming pointlike diquarks. Note that such a decomposition in Dirac 
and Lorentz space holds for the vertex function P) as well. 



The Faddeev solutions are normalized by the canonical condition 

g? 4 p f d 4 p' 



M n A + 



(27T) 4 



(2vr) 4 
d 



(29) 



P "op-G~\p\p,P) 



P=Pn 



3.2 Delta 



The effective multi-spinor for the delta baryon representing the BS wave function can be charac- 
terized as ty*^(p, P)u u (P) where u u {P) is a Rarita-Schwinger spinor. The momenta are defined 
analogously to the nucleon case. As the delta state is flavor symmetric, only the axialvector 
diquark contributes and, accordingly, the corresponding BS equation reads, 

J-^G A 1 (p,k,P)^(k,P) = 0, (30) 

where the inverse quark-diquark propagator in the A-channel is given by 

G a \pXP) = (2n) 4 8"( P -k)S- 1 (p q )(D^)- 1 (p d ) + 

X»'(pl)S T (q)r(pl)- (31) 
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Set I II III IV V 



w 1+ [GeV 2 ] 0.4 0.6 0.8 1.0 1.2 

w 0+ [GeV 2 ] 0.21 0.27 0.32 0.36 0.39 

m l+ [GeV] 0.92 0.91 0.89 0.88 0.87 

m + [GeV] 0.75 0.77 0.80 0.84 0.86 



Table 1: Five parameter sets which describe the physical masses of N, M n =940MeV, and A, 
m A =1230MeV. 

The general decomposition of the corresponding vertex function obtained as in eq. (f23|) 
by truncating the quark and diquark legs of the BS wave function reads 

Ma iMa Ma /Ma 

1 1 Ma iM A iM A K M A %M A 



A^=A+ ( ^ _ _ 7 f 7 " + ___ _ I L_ L) (33) 



Here, is the Rarita-Schwinger projector, 

1 „ „ 2P»P U ipw-pw 

— T T H ! 

3 ' ' 3 Ml 3 M A 
which obeys the constraints 

pn A n» = 7 m A ^ = Q (3 4 ) 

Therefore, the only non-zero components arise from the contraction with the transverse relative 
momentum p 1 ^ = p^ — P^(p ■ P). The invariant functions Di and E{ in eq. (3~2) again depend 
on p 2 and p ■ P. The partial wave decomposition in the rest frame is given in ref. |3( 



3.3 Numerical solutions 



The Faddeev equations for N and A are solved in the baryon rest frame by expanding the 



unknown scalar functions in terms of Chebyshev polynomials of the variable p ■ P |3J] . Thus 
the equations are reduced to a system of homogeneous one-dimensional integral equations. 
Iterating the integral equations yields a certain eigenvalue which by readjusting the parameters 
of the model is tuned to one. As remarked earlier, we are left with one free parameter which is 
taken to be the width W\+ of the axialvector diquark. 

In Tab. [I] we show five parameter sets which lead to a bound nucleon and delta with the correct 
physical masses, M n = 0.94 GeV and Ma = 1.23 GeV. Instead of the A parameters the pole 
locations in the diquark propagators are tabulated which have a more intuitive interpretation. 

We need always a broader axialvector diquark (in momentum space) to fit both nucleon and 
delta. The resulting diquark masses, notably the mass difference m 1 + — m +, agrees approx- 
imately with previous rainbow-ladder results || ^ only for the first two sets with smaller 
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diquark widths. The lattice results of ref. 1 1 1] give m + = 0.83 GeV and m x + = 0.9 GeV, 
within the spread of our parameter sets. 



As mentioned in the introductory chapter, pions are expected to reduce the nucleon mass by 
at least 200 MeV f2~2]| . The corresponding mass shift for the delta will be lower. Therefore we 
investigate the Faddeev equations for nucleon and delta core states of mass 1.2 and 1.4 GeV 
respectively and present the results in Appendix [A[ 



4 Nucleon form factors 



4.1 Electromagnetic gauge invariance 



To calculate form factors, we apply the gauging formalism of ref. [|18| which basically consists 
in coupling the photon to all elements in the kernel G -1 of the nucleon Faddeev equation (|24|). 
Therefore we need the photon vertices with quark, diquarks and the quark exchange kernel. 
Each vertex has to satisfy its WT identity. 



For the quark-photon vertex = r^ BC + r^ T the construction of the longitudinal part, T^ BC , 
which is fixed by the WT identity has been long known [B2). It is given by 

rJ,Bc(*, P) = -vf^±^- - i(p + kr^AA - (p + k)»AB (35) 

where AX = (X& — X p ) / (k 2 — p 2 ) and X& = X(k 2 ), (X = {A, B}). The remaining transverse 
part, is yet undetermined. To ensure multiplicative renormalizibility at the one-loop level, 
an ansatz for this part was proprosed in ref. |33j but it modifies our results for form factors only 



on the level of one per cent. The transverse part might also receive dynamical contributions 
from the p — u meson poles in the q — q vector channel | 3"3|| . In Appendix [FJ, we derive a 



parametrization of such ontributions which is well-constrained by the pion form factor. It is 
given by 



IWM = rf p ^ p e V ^, (36) 
^t-1-69— niQn , (37) 



u p J 0.139 

Q = k-p, q={k + p)/2 < = u" -Q^vQjQ 2 . 

The constants appearing herein are: p mass and decay constant m p = 0.77 GeV and f p = 0.215 
GeV (calculated), a = 0.652 and u 2 = 0.35 GeV 2 . The structure M represents a properly 
normalized vector meson BS amplitude (see Appendix |B|). In eq. ([36]), the combination of the 
exponential and the propagator-like denominator parametrizes the effects of the propagation 



of an off-shell, composite p — u>. For details, see ref. [34 



The diquark-photon vertices and receive contributions from four different diagrams, 
depicted in Fig. H- Besides the photon coupling to the quarks within the loop we need seagull 
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Figure 3: The photon-diquark vertex, here for the scalar diquark. For the axialvector diquark, 
X 5 and M 5m have to be replaced by their corresponding counterparts. Note that in actual 
numerical calculations the first two diagrams are equivalent if the quarks are identical. 



graphs which describe the photon coupling to the vertices x an d X M to recover the WT identity 
(k — p) M Fg + j 1+ j = (D^'Y 1 (k) — (.D^) -1 ^). The functional form of these seagull vertices has 
been derived in ref. ]19| and they read 

(4// - QY 
4 P '-Q-Q 2 



(M a r(p',Q;q a ,qa) = q a [xV - Q/2) - X a (p')} + 



% W-Q + Q 2 [xa{p ' + Q/2) " xV)] • (38) 

The photon momentum is denoted by Q = k — p. The relative momentum p' between the two 
quarks with charges q a and qp has been defined in eq. @. The conjugated vertex is obtained 
by replacing x - * X, Q — Q an d interchanging q a <-> g^. 

Photon-mediated transitions between scalar and axialvector diquarks are also possible. The 
corresponding (anomalous) vertices r +_x+ resp. rx+_ + are described by diagrams like the 
first two in Fig. [| Seagulls give no contributions to these vertices. 

The quark exchange kernel, as given in eq. (0), consists of expressions 

~ X a S T x b ( a ; b = {5, u}). Complete gauging leads to a diagram where the photon couples to 
the exchange quark and two diagrams where the photon couples to the vertices x a an d X b - The 
latter couplings are described by the vertex as given in eq. d38|). The proof that the gauged 
quark exchange kernel obeys its WT identity can be found in ref. [Z9|. 



In summary, to obtain the complete nucleon current matrix element, we have to calculate the 
diagrams shown in Fig. |j. We remark in passing that the normalization condition for the 
nucleon Faddeev amplitudes, eq. (|29|), is only compatible with the correct nucleon charges, i.e. 
G\e(0) = 1 (proton) and G\e(0) = (neutron), if all diagrams of Fig. § are taken into account. 



12 





Figure 4: Nucleon current matrix elements: (a) Impulse approximation diagrams. The photon- 
diquark vertex consists of the elements given in Fig. |3|. (b) Exchange kernel diagrams. 

4.2 Numerical calculations 

We extract the Sachs electromagnetic form factors from the current matrix elements by the 
following traces: 

G e(Q 2 ) = ^sT±{P f \J''\P i )pP, (39) 

i M 2 / * \ 

G M (Q 2 ) = -^Tr(P f \r\P t )(Yh, (V)t = 7 m -^) • (40) 

Here, P = (Pi + -P/)/2, and the initial and final states |Pj) and (Pf \ are given by the numerical 
solutions for the matrix valued wave functions ^(p,Pt) and ^(k,Pf), cf. eqs. (|2^ , |27|j28|) . 

Due to the complicated singularity structure of the single diagrams and due to limited computer 
resources we obtained fairly accurate numbers for the form factors only up to momentum 
transfers of Q 2 = 2 GeV 2 . These technical obstructions do not interfere with the conclusions 
we will draw, though. A detailed discussion of the technicalities is deferred to Appendix 0. 

We have calculated the form factors for the five sets tabulated in Tab. The results for the 
electric form factors of proton and neutron are shown in Fig. [5[ The proton Ge (left panel) 
becomes steeper with decreasing diquark widths W\+ and Wq+. This is in agreement with 
intuition since the scalar diquark correlations give the most important contributions to this 
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Figure 5: Depicted are the electric form factors of proton and neutron. Experimental data for 
the proton are taken from the analysis in ref. [Q. For the neutron, data are from ref. p6 
(diamond), from ref. |37| (square), from ref. ]38| (triangle up) and finally from ref. |39| (triangle 
down) . 



form factor, and with decreasing width also mn+ decreases, so that these correlations become 
wider in the "center-of-mass" position variable as in the relative position variable. However, 
we also observe for Set I an interesting deviation of the form factor from the dipole shape which 
will be discussed below. 



Let us turn to the results for Ge of the neutron (right panel in Fig. ||). All data sets predict a 
positive form factor which is slowly falling for larger Q 2 > 1 GeV 2 . No data set can reproduce 
the experimental neutron charge radius or come close to it. This is in remarkable contrast to the 



results in refs. [15, 16] where rather simple parametrizations of the q — q t matrix were employed. 
As for ref. JT5| , the good description of Ge was mainly a result of the cancellation between 
quark and diquark impulse approximation diagram. The former contribute negatively, the latter 
positively and by virtue of the simple approximation of the q — qt matrix by free spin-0/spin-l 
particles and the corresponding free photon vertices, the diquark contributions fall slower and 
thus render Ge positive. In this study, we have resolved the diquarks (see Fig. (d)) and the effect 
of diquark contributions falling more slowly is almost absent, thus the charge radius becomes 
very small. Only with a proper resolution of the diquarks an asymptotically correct description 
for the form factors is possible at all, thus we conclude that the neutron charge radius must 
be accounted for by other mechanisms such as a neutron dressing by pions. Nevertheless the 
positivity of Ge for higher momentum transfers is a result of the fully relativistic treatment. 



Set 


I 


II 


III 


IV 


V 


fi p [n.m.] 
H n [n.m.] 


3.05 
-1.78 


2.94 
-1.65 


2.86 
-1.55 


2.79 
-1.47 


2.73 
-1.40 



Table 2: The magnetic moments of proton (fi p ) and neutron (/i„) for the five parameter sets 
given in table 1. 
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The Ratio |iG E / G M (Proton) 




Figure 6: Magnetic form factors and Ge/Gm- Experimental data for proton's Gm are from 
ref. |35) and data for Ge/Gm have been reported in ref. (40] . 



Fig. || left panel, shows the nucleon magnetic form factors and Tab. |2] the corresponding 
magnetic moments. Also the magnetic radii become larger with decreasing diquark widths, as 
well as {ip and \n n / n P \. 

Our results for the ratio ^pGe/Gm, currently under intensive experimental scrutiny, are shown 
in the right panel of Fig. |j. Although the results consistently put that ratio below 1, the 
available experimental data are underestimated considerably. Even worse, as going towards 
more realistic electromagnetic radii (with decreasing set number) the ratio becomes smaller 
and smaller. Before giving a reason for the underestimation of [i p Ge/Gm, we will examine the 
influence of the vector mesons in the quark-photon vertex. 

Full vertex vs. Ball-Chiu vertex 

Since the resonance contribution is ~ Q 2 near Q 2 = 0, it does not give any contributions to 
the magnetic moments, thus these are accounted for by the Ball-Chiu vertex alone. As for the 
pion form factor, to which the resonant vertex was fitted, it does give sizeable contribution to 
the charge radii. It amounts to 21-23 % in the case of {r p )\ for all data sets and is therefore 
of the same relative size as the contribution of the resonance to the pion charge radius, see 
Appendix ||. The neutron charge radius is usually a bit smaller, as the electric form factor is 
quenched a bit more when the full vertex is employed, nevertheless the differences are small 
and can not account for the discrepancy with the data. 

In Fig. [j] we show the influence of the resonance term on proton's Ge and the ratio Ge/Gm 
for set II. We investigated two cases of resonance contributions, the full vertex which includes 
the transverse terms as in eq. [36] and a transverse vertex which includes only the leading 7^ 
term. In both cases, the vector meson amplitude is normalized and the damping constant a 
has been fitted to the pion form factor. For small Q 2 , both resonance parametrizations lead 
to a Ge falling more quickly, with almost no quantitative difference. For intermediate Q 2 , the 
subleading vector meson amplitude ~ gx leads to a slightly enhanced Ge- This effect is visible 
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Figure 7: Comparison between results for the full vertex, the Ball-Chiu (BC) vertex plus a 
dominant resonance term, and the pure Ball-Chiu vertex. In the first two cases, the resonance 
terms have been fitted to the pion form factor. 

more clearly in the ratio Ge/Gm- The subleading amplitude also quenches Gm a bit such 
that the ratio comes out more than 10 % larger than for the BC vertex in the intermediate Q 2 
domain. 

We see that the vector meson resonance has a sizeable influence on the proton charge radius 
(of the order 1/4), and subdominant amplitudes of the vector meson can influence Ge/Gm by 
10-15 %. Nevertheless for this observable a discrepancy remains, and a reason can be found 
by analyzing the Ball-Chiu quark-photon vertex. This vertex always appears with quark legs 
attached, i.e. in the combination = S(k)r^ BC S{p). We rewrite in the following way: 

r{? = ff,BC + It ■ ( 41 ) 

The vertex with legs fulfills the WT identity Q^T^ = S(p) — S(k) (note that the propagator 
and not its inverse appears on the r.h.s.). The term Fg BC is constructed via the Ball-Chiu 
technique to satisfy this identity: 

f?,Bc(*> p) = -ir avk ^ avp - up + kr^±l Aav + {p+ krAas . (42 ) 

The remainder, T^, is transversal, and after some Dirac algebra one finds 

f£ = [(T-^g+iAM^-^)-Q(k+ P r)]^^^- 

G -f%T^ k ^ ~ Qr + iAM ^ -P 2 )-$( k + PT)} • (43) 

Due to the running mass function, the terms proportional to AM = (M(k 2 ) — M(p 2 ))/ (k 2 —p 2 ) 
are non-zero. Precisely these terms give a fairly large negative contribution to the proton's 
G E: thus causing a deviation from the dipole shape. This effect is absent for G M - Of course, a 
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Figure 8: Left panel: lattice data from ref. |^5[ for the quark mass function, extrapolated to 
the chiral limit. Right panel: the ratio Ge/Gm for calculations with modified transverse part 
of the quark-photon vertex. 



dynamic quark mass function of the kind depicted in Fig. [I] is phenomenologically required and 
thus, as in the case of the neutron electric form factor, the effect of other contributions besides 
the quark valence core should be sizeable for this observable. Nevertheless, some parts of the 
discrepancy (though not all) could be attributed to a stronger influence of the sub dominant 
vector meson amplitudes. This amounts to a shift of the off-shell vector meson contributions 
between the Ball-Chiu vertex and the transverse contributions and thus between the quark 
propagator and the off-shell vector meson amplitudes. 



To get a first estimate how a change in the quark mass function might influence the ratio 
Ge/Gm, we recalculated the form factors using the wave functions from parameter set II, but 
replaced all occurences of AM in the expression for Hf, eq. fl43|), by (a) zero (corresponding to 
a momentum-independent quark mass), and (b) by a fit to the most recent lattice data [[45 



Since we are using wave functions calculated with the quark propagator parametrization from 
eqs. (0,|), this procedure is somewhat inconsistent, but may give qualitative indications to the 
behavior of Ge/Gm- Nevertheless, gauge invariance remains intact since the AM terms only 
appear in the transversal part of the vertex. We remarked earlier that the quenched-QCD 
(Landau gauge) lattice data seem to suggest a somewhat broader quark mass function than 
the parametrization employed here. This can be seen from Fig. |] (left panel) where chiral limit 



extrapolations of the data and a fit from ref. [45| are given. Since the functional form used in 



the lattice fit, M(p 2 ) = cA 1+2a /(p 2a + A 2a ) with a = 1.52, has a cut along the negative p 2 axis, 
it is not well suited for our calculations. We chose to re-fit the data to the entire function 

M(q 2 ) = Ci p-exp(-g 2 2 -^) \ 2 ^ [?2=p2/(1GeV2)] _ (44) 

For the choice of parameters C\ = 0.4 GeV and c 2 . = 0.45 our fit is also depicted in Fig. [| 
We remark here that both fits, lattice and entire one, lead to a 20 % underestimation of the 
pion decay constant f„ (using the formula from ref. |17|]), reflecting the uncertainties in lattice 



extrapolations and the shortcomings of the quenched approximation. 
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Results for the ratio Ge/Gm are plotted in Fig. || (right panel). The ratio is lowest for the 
consistent calculation. If eq. (fPJ) is used in the vertex T^, the curve is shifted upwards and 
the proton magnetic moment is somewhat smaller. Surprisingly the simple constituent quark 
assumption AM = delivers the best results compared to the experimental data, though 
employing the approximate wave functions prevents us from drawing precise quantitative con- 
clusions. Nevertheless we see that the ratio Ge/Gm is sensitive to the precise form of the 
running quark mass. Furthermore we note that nearly the whole effect comes from the quark 
impulse approximation diagram (the second one in Fig. ^). 

Clearly more precise QCD lattice data for the quark propagator and/or DSE/BSE studies of 
the three-quark and quark-photon systems are desirable. 

Form factor results for nucleons with higher quark-diquark core mass are presented in Ap- 
pendix |A[ There we find that the neutron electric form factor does not change, and that the 
even steeper mass function in the time-like domain deteriorates the ratio Ge/Gm even further. 



5 Summary and conclusions 

In a step towards the solution of the full covariant Faddeev equations for baryons, we have mod- 
eled two-quark correlations by assuming them separable and by summing quark polarization 
loop diagrams. In this case, the Faddeev equations reduce to a Bethe-Salpeter equation which 
has been solved exactly. The technique employed in calculating the two-quark correlations 
effectively reduced the number of model parameters to one, the diquark width. 

The nucleon form factors have been calculated in a scheme which preserves the Ward-Takahashi 
identities for the basic two-point function, the quark propagator, for the four-point function, 
the quark-quark scattering kernel and finally for the quark-diquark kernel of the Faddeev 
equations. Consequently the current is conserved. Constrained by the pion form factor, effects 
from vector mesons have been included in the quark-photon vertex. 

Results reveal two effects. If the proton electric and magnetic radius is to be described correctly, 
the ratio Ge/Gm is severely underestimated. This is a consequence of the parametrization of 
the dynamic mass function of quarks in accordance with results from Dyson-Schwinger and 
lattice studies. Furthermore the substructure of the two-quark correlations which is resolved 
by the photon renders the neutron form factor positive, but quite consistently for all parameter 
sets the corresponding charge radius cannot be described. Assuming core nucleon states with 
higher mass does not alter the above findings. 

Certainly the precise shape of the form factors is expected to vary if the technical simplifi- 
cations can be rendered obsolete, such as the separability of the two-quark correlations and 
the treatment of the vector meson contributions to the quark-photon vertex. Nevertheless it 
seems possible that the qualitative features will remain valid, i.e. the quenched neutron electric 
form factor for a correctly resolved q — q matrix and the underestimation of Ge/Gm due to 
the quark-photon vertex with running mass function. As described, the vector meson contri- 
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butions cannot compensate this effect. The running quark mass and vector mesons are usually 
not considered in non- or semirelativistic quark models and urge us to a cautious interpretation 
of corresponding results, see e.g. ref. |4l| . Thus the investigation presented in this paper point 
towards the necessity to incorporate non-valence quark physics into the description. Covariant 
studies of the effect of e.g. the pion cloud within covariant bound-state perturbation theory 
are clearly desirable. 
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A Results for core nucleon and delta states with higher 
mass 



Pion cloud corrections will lower the mass of the nucleon by more than 200 MeV. The corre- 
sponding mass shift of the delta will be somewhat lower as the nucleon-delta mass difference is 
partly also a consequence of pionic dressing. For exploring effects of this shift on observables, 
we fix the core masses to M n c = 1.2 GeV and M^ c = 1.4 GeV. Possible parameters which lead 
to a solution of the Faddeev equations are given in Tab. [| We note that less than half of the 
mass difference between core and physical nucleon state can be attributed to the shift in the 
diquark masses which are about 100 MeV larger than those in Tab. [IJ If one assumed pointlike 
diquarks, the scalar diquark would not be changed by pionic corrections at all, but the quark 



substructure allows for pion dressing [42 



The proton magnetic moments are slightly reduced compared to the results for nucleons with 



Set I II HI IV V 

wi+ [GeV 2 ] 0.4 0.6 0.8 1.0 1.2 

w + [GeV 2 ] 0.21 0.27 0.32 0.37 0.41 

m 1+ [GeV] 1.02 1.01 1.01 1.00 0.99 

m + [GeV] 0.89 0.90 0.91 0.92 0.94 



Table 3: Five parameter sets which describe core masses of nucleon and delta, M n c = 1.2 GeV 
and M AjC = 1.4 GeV. 
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M n [GeV] Set 


I 


II 


III 


IV 




[n.m.] 0.94 


3.05 


2.94 


2.86 


2.79 




1.2 


2.88 


2.80 


2.73 


2.69 




[n.m.] 0.94 


-1.78 


-1.65 


-1.55 


-1.47 




1.2 


-1.79 


-1.69 


-1.61 


-1.56 



Table 4: Nucleoli magnetic moments, compared betweens sets with physical nucleon mass and 
core mass. For each pair of these sets, W\+ is equal and w + differs by at most 5 %. 




0.5 1 1.5 2 0.5 1 1.5 2 

Q 2 [GsV ! ] Q* [GeV ! ] 



Figure 9: Form factor results compared for data sets describing nucleons with physical mass 
(full line) and core mass 1.2 Gev (dashed line). Left panel: proton's Ge and [iGe/Gm- Right 
panel: neutron's Ge- 



physical mass whereas the neutron magnetic moments come out larger by a few percent, see 
Tab. £|. We present exemplary results for the form factors Ge and proton's Ge/Gm in Fig. || 
comparing Set III from Tabs. [I] and |3|. The characteristics found here apply to the other 
data sets as well. The proton electric form factor becomes steeper thereby approaching the 
experimental curve and at the same time the ratio Ge/Gm moves even further away from the 
data. Since for higher core mass the quark-photon vertex tests the running mass function 
deeper in the timelike domain, its steeper derivative there leads to this effect. The neutron 
electric form factor remains unchanged. 



B Resonance contribution to the quark— photon vertex 

In this section we describe shortly the procedure to fix a p — uj resonance term in the quark- 
photon vertex. A somewhat longer discussion of the subject and the techniques used herein 
can be found in ref. [TJ. 

The transverse part of the quark-photon vertex, T will certainly receive resonance contribu- 
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tions from the p — lo mesons. Assuming isospin symmetry and neglecting the decay width, the 
full vertex is dominated near the resonance by the term 

q Q 2 + m 2 

Q = k-p q={k + p)/2. (46) 

Here, m p = 0.77 GeV and f p = 0.216 GeV are the mass and electromagnetic decay constant of 
the p meson. The BS amplitude M (g; Q) of the p meson is transversal (Q ■ (ft = 0) and obeys 
the canonical normalization condition 

2Q^ = I 3 Tr J 0^ 4 nq; Q) S(-Q/2 + q) Q) ^ ( ^ + g) , (47) 

if we assume that an interaction kernel for the corresponding BS equation is independent of Q 
(as e.g. a dressed gluon exchange between the quarks). The factor 1/3 comes from the sum 
over the three p polarizations, and the factor 3 is the result of the combined flavor and color 
trace. In color space, both p and u amplitudes are the unit matrix 5ab, i n flavor space we have 
(r 3 ) ab /V2 (p) and (r°) ab /V2 (u). 

Away from the resonance mass-shell the corresponding contribution to the quark-photon vertex 



is not fixed uniquely. The most thorough study of it is ref. [|34|], which calculates the (gluon- 
ladder) dressed quark-photon vertex for the evaluation of the pion's form factor. The findings 
of ref. |34| may be neatly summarized by the following points: 



• ^(Q'iQ) ~ ^7t ^U 1 ? 2 ) + ^5(<? 2 ) represents a good approximation to the BS solution 
for the p meson ( vt^ = v M — Q^v ■ Q/Q 2 )- It reproduces the mass and decay width within 
5%. 

• The Dirac structure accounts also for the bulk of the resonance contribution to the 
pion form factor, whereas terms ~ Qt provide corrections to these contributions on the 
level of 10 %. Using the p BS amplitude off its mass shell in the manner of eq. fl45|) gave 
a good approximation to the quark-photon vertex resonance contributions. 

• The Dirac structure (qt — 7t^t)^* becomes more important for intermediate Q 2 in the pion 
form factor but the off shell extrapolation of the corresponding BS amplitude structure 
proved to be difficult. Thus we neglect this term. 



We therefore adopt an off-shell parametrization of the resonance term in the quark-photon 
vertex, 

T£t = ni) "7^7t2~7~ — 2 • (48) 

Near the p mass shell, Q 2 = —m 2 , eq. fl48D reduces to eq. (f45|). The exponential ensures that 
for high spacelike Q 2 the resonance term vanishes and the quark-photon vertex reduces to the 
Ball-Chiu vertex, eq. (|35|) . The resonance term vanishes also for Q = as it should be, since at 
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this kinematical point the quark-photon vertex is completely fixed by the (differential) Ward 
identity 



We model the p BS amplitude close to the results of refs. |34|, [|3|] by employing the one- 
parameter form 



(49) 



The normalization constant N p is implicitly given by eq. fl47|) and the width parameter uo p will 
be fixed by the experimental value for the decay constant, whose theoretical expression is 



f P = — Tr 

m„ 



d A q 



-i^) S(-Q/2 + q) <ff{q) S(Q/2 + q) 



(50) 



=(Q,im p ) 



The only unknown parameter which remains in eq. (fiSp is the constant a which describes the 
damping of the off-shell resonance contribution. We fit it to the pion form factor in the range 
Q 2 = [0, 1.6] GeV 2 where the experimental data are well described by the monopole fit [55 



0.529 GeV^ 

In impulse approximation, the pion form factor is given by 

Hi 



FAQ 2 ) = K^— ■ ( 51 ) 

1 + 



UQ 2 ) = ^Tr J^Mk^Si^M^Sik^T^k^Sik^ (52) 



k + = P/2 + k + Q/2, k f = k + Q/A, 
k_ = p/2 +k- Q/2 , k t = k- Q/A , 
ki = -P/2 + k, 



Q = (0,0,|Q|,0), P= (0,V< + Q 2 /4) • 

For the region of momentum transfer in consideration, the truncation of the pion's BS amplitude 
to the leading amplitude which is determined by chiral symmetry, 

MP 2 ) = ?P- 75 , (53) 

is an excellent approximation. Indeed, the quark propagator used herein has been fitted to just 
give N n = f n = 93 MeV, as expected physically^. 

As a result, we obtain p BS amplitude width of w 2 p = 0.35 GeV 2 which in turn gives f p = 0.214 
GeV. The damping factor a = 0.652 results in a pion form factor as shown in Fig. [10], compared 
to the result obtained with only the Ball-Chiu vertex. We remark that the latter yields a pion 
charge radius of r 2 = 0.31 fm 2 , whereas with the resonance contribution we obtain r 2 = 0.44 
fm 2 , in acordance with the experimental value. Thus, about 30 % of the charge radius is 

1 We normalize the pion BS amplitude in color space by 8ab and in flavor space by (T k ) a b 
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Pion Form Factor 




Figure 10: The pion form factor, calculated with the Ball-Chiu vertex (thin line) and the 



vertex including the vector meson resonance (thick line). Experimental data are from ref. p6 
(diamonds) and ref. |Q (filled circles). 



attributed to the vector mesons. This contribution is only half of the value obtained in ref. 33 



and thus reflects the model dependence of the off-shell extrapolation: The parametrization 
for the quark propagator used here has a steeper Z(p 2 ) = 1/A(p 2 ) than the corresponding 
renormalization function of the quark propagator obtained in ref. The Q 2 — variation of 

the scalar function multiplying the dominant Dirac structure ~ 7^ in the quark-photon vertex 
is therefore already rather steep for the BC vertex used in this work and the resonance part is 
seen too influence the form factor only up to Q 2 = 1 GeV 2 . 



C Singularity structure of the form factor diagrams 



Singularities in the diagrams are present through the quark-photon vertex T^, see eq. (|35|). It 
contains the scalar functions A and B which are defined by 

A ^ = MP 2 ) P 2 + M 2 (p 2 ) ' (54) 
B(?2) = S^iW ' (M(p 2 )= ff . s( P 2 )KV)). (55) 

We see that these functions have poles whenever ay or p 2 + M 2 have zeros. The poles in 1/ ay- 
do not matter since in the current matrix element diagrams, always appears with quark 
legs, ST^S, and these legs cancel the pole. Such a mechanism is not present for the poles in 
l/(p 2 + M 2 ). A numerical search revealed the following poles, being closest to the origin in the 
complex p 2 plane: 
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Re p 2 Im p 2 

GeV 2 
-0.067 ± 0.207 
-0.167 ± 1.116 
-0.224 ± 1.360 
-0.293 ± 0.742 



We see that the pole locations appear in complex conjugate pairs. It is only the first pair of 
poles in the list which will have an impact on our calculations. To see that, let us consider the 
single form factor diagrams. 

We calculate the form factors in the standard Breit frame where 



P = (0,0,0,V^ 2 + Q74) 



Pi = (0,0,-|<2l/W^ 2 + <274) (56) 
P f = (0, 0,+\Q\/2,i^M 2 + Q 2 /A) 

We start with the current matrix elements of the impulse approximation (diagrams (a) in 
Fig. |J). The quark diagram (where the diquark is spectator) is given by 

Pi = k-(l- i])Q/2 , k q = r]P + k + Q/2 , 
Pf = k + {l- i])Q/2 , Pq = v P + k-Q/2 , 

k d = (l-T])P-k. 

The loop momentum k is real, but the quark momenta k q and p q at the vertex are not. Since 
their imaginary part, r]P grows with increasing Q 2 , the integration domain will cross the pole 
locations p 2 ole . The limit for Q 2 such that the integration domain will be free of these poles is 

g 2<2 bpolc|-RePpole _ 4M 2_ (5g) 

Beyond this limit, the integration path in the variable k± has to circumvent the poles or, alterna- 
tively, the sum of a principal value integral with the original path and a closed contour integral 



encircling the singularities must be calculated. This fact has been overlooked in ref. [16]. While 



this procedure has been carried out for real poles in ref. |T{|, for complex poles the knowledge 



of the wave function at complex relative momenta between quark and diquark is required. To 
obtain the wave function at these points is in principle possible, but the implementation in a 
current matrix element code is at present not feasible. 



Choosing small momentum partitioning parameters rj shifts the limit (|58f) to larger values. 
Due to the presence of the diquark pole, solutions of the Faddeev equation are restricted by 
rj > 1 — m +/M n , in practice we have to restrict ourselves to r] = 0.32 to have the Chebyshev 
expansion of the wave function converge for both the Faddeev solution and the calculation of 
the current matrix elements. This yields the limit Q 2 < 2 GeV 2 , by virtue of eq. (]58f). 
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In the diquark diagram (with the quark being spectator), singularities occur only where the 
photon resolves the quark loop, i.e. in the integrals 



W = f^^\Vf,Pf)S-\k q )Y^{k dlPd )^\ Pl) P t ) (59) 

Pi = k + 7]Q/2, h d = (l-ri)P-k + Q/2, 

Pf = k- V Q/2, Pd =(l- V )P-k-Q/2, 

k q = f]P + k 

T» + (k d ,p d ) = Tr l^f( q+Q/ ,)S(q 2 )T^q 2 ,q 1 )S(q l ) X 5 ( q -Q/^S T (q 3 ) 

(60) 

9i = (p d + k d )/A + q-Q/2 , q 2 = [p d + k d )/A + q + Q/2 , 

<?3 = (Pd + k d )/A- q . 

While these equations describe the contribution of the scalar diquark, similar expressions hold 
for the axialvector diquark and the scalar-axialvector transitions. The imaginary part of the 
quark momenta q\ and q 2 which enter the quark-photon vertex is given by (l — T))P/2, therefore 
we can apply eq. ( j58|) for a pole-free integration domain upon the replacement 77 — ^ (1 — rj)/2. 
If we want to calculate pole-free up to Q 2 = 2 GeV 2 , we find r\ > 0.36. Thus the Faddeev 
solutions and these diagrams have to be calculated with a different momentum partitioning 
parameter than the Faddeev solutions for the quark diagram. This is of course possible since 
the Faddeev solutions have been obtained fully covariantly, i.e. the full dependence of the wave 
function on p 2 and p ■ P has been retained (p is the relative quark-diquark momentum and P 
is the total nucleon momentum). 

Singularities in the exchange kernel contributions (diagrams (b) in Fig. f|) are present in the 
diagram where the photon couples to the exchange quark. The corresponding current matrix 
element is given by 

(JQ = /^/^^Wfe)(%i)^(gi,g 2 )% 2 )) T x 

X 6 (fc.)* b (p,Pi) (a, 6= {5,fx}) (61) 
qi = - p -k+(l-2r))P-Q, Pa = k + Z-{l-3ri)^-(7i-l)Q 

q 2 = - p -k + (l-2ri)P + Q, fc s = p+^-(l-377)^ + (77-l)| 

Here the imaginary part of the momenta q\ and q 2 (appearing in the quark-photon vertex) is 
given by (1 — 2r])P. Using eq. (|58|) again (with 77 — > (1 — 2rj)), we find the condition 77 > 0.34 
if we want to calculate this diagram without encountering poles up to Q 2 = 2 GeV 2 . 
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